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Abstract 

Why did mankind create, and continues to create, societies? What is it that seems to keep 
societies together? And then, are there natural boundaries for societies mankind would not 

exceed? 

The first question is rather philosophical, and we comment only briefly on it. The second 
one displays several quantifiability aspects which would be open to statistical hypothesis testing, 
but we only treat it in as much as it concerns the third question, which attracts our particular 
interest. Here we shall provide a mathematical answer obtained from what wc propose as a 
global mathematical model for societies. It is built on branching processes, and submitted to 
two natural hypotheses. Still rudimentary, our model allows nevertheless to take into account 
essential features of life within a society. Apart from reproduction of individuals, it incorporates 
the following factors and interdependencies: the desire to survive, heritage and production of 
resources, consumption of resources, policies to distribute resources among individuals, and, 
decisively, the right of emigration. We look at different submodels of the model, characterising 
different societies. These are defined by the type of control they exercise through different 
policies to distribute resources among their individuals. 

Our main results are: Firstly, in the global model we consider, the answer to the third 
question is affirmative. Secondly, we can determine these boundaries. Thirdly, we can also 
determine their critical parameters, and interpret them. Indeed, there are exactly two societies 
which form an envelope in the sense that all societies have to live in the long run between 
these two boundaries. We call the boimdaries weakest-first society and strongest-first society. 
Interestingly, each is related with a society form we believe to recognise rather well. It seems 
that, with respect to societies, mankind has already come close to testing the limits. 
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1 Introduction 

Different theories in anthropology and sociology suggest that man has a natural tendency for social- 
ising. Some of these theories see this tendency as a fact which is, a priori, independent of the idea 
of synergy effects and utility. Others contrast in the sense that they consider the essential reasons 
why people would join together to form groups or larger communities lie in the desire to profit from 
an efficient collaboration, as e.g. in hunting, defense, agricultural activities, industrial production, 
etc. As so often, the truth may lie somewhere in the middle. But then for any theory explaining 
the phenomenon of socialising, other questions come up naturally: Why do groups and larger com- 
munities have a tendency to go further and develop "structured" societies with, often enough, very 
restrictive rules? And then, how far can man go? Are there natural boundaries for societies? 
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Certain aspects of these questions are equally interesting for animal societies. In fact, throughout 
this paper we shall always speak of "individuals" to make clear that, although the motivation stems 
from thinking about man, wc keep general populations in mind. 

The present paper tries to add to our understanding by looking at certain questions which nat- 
urally arise from basic mathematical models. We link a population model (in terms of a branching 
process model) with different rules to distribute resources without which the members of the popula- 
tion cannot live and reproduce. The philosophy behind the model we propose is that if mankind has 
a natural tendency for socialising, then it must also have the desire to see this social environment 
survive;. With limited life time for each individual, survival means survival of their descendants. 
Moreover, to maintain a society, sufficiently many individuals in the population should have access 
to a minimum standard of living. 

Theoretical desiderata must here be adapted to reality. So, for instance, in branching process 
models, survival is usually seen as survival "forever" . We know that with our current understanding 
of life in our solar system the "forever" is bound to have the meaning of "very long." Also, in 
population models based on branching processes, a typical feature is that the population will either 
get extinct or, alternatively, explode, the latter being again incompatible with reality. Nevertheless, 
similarly as in modelling physical phenomena, where extrapolation often makes no sense outside 
certain thresholds but gives interesting insight within certain bounds, our models may be seen as 
adequate as a first approach. 

The first-named author has been thinking about such problems repeatedly for many years. Much 

effort has gone into constructing a tractable model. He had spoken about resource dependent branch- 
ing processes as early as 1983, given a second talk on this around 1995, and a third one around 2002. 
Seeing his results as being immature, he had never submitted any of them for publication. As far 
as we are aware, although in the meantime several interesting new branching process models have 
been created, this type of interdependent models has not been studied. Hence the models and results 
presented in this paper are believed to be new, and, now refined, of true interest. 

2 Objective and content of this paper 

Our objective is to answer the questions we asked through a suitable modeling and mathematical 

conclusions. Our approach is to construct branching process models in which individuals have to 
create resources in order to be able to live and where the control is implemented by policies to 
distribute resources among individuals. 

The distribution of resources is effected in each generation, the planning is for the longer-time 
future development of the population, however. We suppose that individuals inherit resources from 
preceding generations, consume resources, and create new resources. The resources an individual 
can use during its lifetime are supposed to determine its standard of living. However, for simplicity, 
we model this relationship between individually available resources and the individual perception of 
standard of living only in a "binary" way, as will be explained below after the main hypotheses. 

The inherited resources, plus the newly created ones, are considered as being the individual's 
contribution to society. We model heritage and production (combined) of resources of individuals as 
independent random variables with a common continuous distribution function. 

Our interest focusses on the implications of the following two basic hypotheses on the policy by 
which the population distributes its resources: 

Hypothesis 1: Populations want to survive; 

Hypothesis 2: Individuals prefer a higher standard of living to a lower one. 

The desire to survive is here understood as the objective to have for the society as a whole a positive 
probability of surviving forever (in the strict sense) . We suppose that populations determine rules 
for distributing resources to individuals who have random resource claims. If certain rules allow for 
a positive probability of survival whereas other rules would not achieve this, then the objective to 
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survive is taking priority, and the rules would be changed accordingly. 

Apart from the reproduction mean per individual and the resource production mean per individ- 
ual, the distribution of resource claims and the rules to try to satisfy these claims will play an equally 
important rule for survival. We think of each claim as being the outcome of what an individual, with 
its own power of conviction, will be able to defend within the society. The random claims of individ- 
uals are also modelled as independent random variables governed by a known continuous distribution 
function. The idea is of course that the latter may be chosen as a function of the distribution of 
resource production. For our objective in the present paper it will however suffice to consider both 
distributions as fixed. 

The effect of the interplay of claims and available resources is modelled as follows. Each individ- 
ual is supposed to emigrate if and only if its individual resource claim is not completely satisfied; 
otherwise it remains a member of the population until the end of the generation. Emigration is 
supposed to happen before an individual produces offspring. Hence the resource assignment (of- 
fered standard of living) is seen by an individual as being either sufficient, implying "stay", or else 
insufficient, implying "leave". This is the "binary" measure of satisfaction we announced above. 

The reproduction of the population itself is assumed to be in accordance to the rules of a Galton- 
Watson branching process (GWP), i.e. members of the population reproduce independently of each 
other with the same distribution of the number of offspring. 

Related work 

Our model is an asexual controlled branching process (BP), where controlled should be understood 
in a larger sense. The general control is governed by a population-size dependent function of sums 
of dependent variables, and self-imposed. This strong dependence property excludes the powerful 
generating function machinery, of course. Moreover, although still rudimentary, the model seems no 
longer tractable for other strong tools as e.g. martingale arguments. 

As we understand, and explain below, our model does not fit other models studied in the literature. 
Nevertheless, related work is sincerely acknowledged for having helped to get a feeling what result 
one can, or cannot, possibly hope for. 

The model is neither a BP with varying environment (see e.g. Cohn [1996]) nor a BP with random 
environment. See Jagcrs [1975] for a clear analysis of the connection between these two types, and 
e.g. Haccou, Jagers, and Vatutin [2007] for newer developments. Our model is no multi-type BP 
model, and neither a pure population size-dependent model. It is a Markov process, as we shall see, 
but no phase- type Markov model or decomposable BP (see Hautphcnne [2012]) can play the control 
we have in mind. 

Early work on controlled branching processes confined interest to control through bounds imposed 
on the growth of GW-type processes. Sevast'janov and Zubkov [1974], Schuh [1976] and others mod- 
ified the number of individuals which are allowed to reproduce in each generation by corresponding 
deterministic functions. Bruss [1978] considered a GWP with a non-specified absorbing process for 
which only the expected influence is known. Yanev [1976] studied so-called ^-branching processes 
where the growth of the GWP reproduction is controlled by random numbers of offspring which are 
allowed to reproduce. A similar model for random control functions was studied in Bruss [1980], and 
in more generality by Gonzalez, Molina, and Del Puerto [2002]. The same authors also examined 
L2-convergence for such processes (see Gonzalez, Molina, and del Puerto [2005]). Population-size 
dependence is another interesting access to control in branching process models. These were studied 
by Klebaner [1985] and Cohn and Klebaner [1986]. Xu and Mannor [2012] proposed a special class 
of controlled branching processes involving the notion of "resources". Motivated by applications 
in marketing, their objective is to control independent subpopulations (multi-type model) in such 
a way that they grow as quickly as possible under an (optimal) assignment of resources. Relative 
frequencies of types were studied in Yakovlev and Yanev [2009]. 
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3 Mathematical Model 



We start by a formal definition of the global model which will define the common structure of all 
the models we consider. Here we postpone the explanations which will become clear as soon as we 
introduce those specific models which attract our particular interest. 

Let (I?^, X^, i?^)„£N fegNg be a double array of triplets of random variables defined on a proba- 
bility space {Q,J^,P), where all triplets (Z)^,X^,i?^) are i.i.d. with the same distribution as some 
triplet of random variables {D,X,R). 

-D^, X*' and i?,'^, ^ ^o) represent the number of offspring, the resource claims and the production 
of resources (respectively) of each individual (labeled by k) in generation n. Here, resource creation is 
supposed to summarise production of resources and non-consumption of resources originally claimed. 

We make several assumptions, which are natural with respect to our interpretation: 

(Al) D is an integer-valued non-negative random variable with law 

Pk = P[D = k], 

satisfying pq > and, to avoid trivialities, satisfying also Po +Pi < 1- Let m denote the mean 
reproduction rate, i.e. the expected number of offspring of one individual; we suppose 

oo 

E[£)] = ^ kpk =: m < oo. 

k=l 

(A2) i? is a real-valued non-negative random variable. Let r denote the average resource creation; 
we suppose 

< E[R] =: r < oo. 

(A3) X is an absolutely continuous real-valued non-negative random variable. Let F be the dis- 
tribution function; we assume that the density / = dF/dx has bounded support [a,b]. Let fj, 
denote the average resource claim of one individual. Hence, 

< E[X] =: iKoo. 

(A4) We shall assume throughout this paper that the random variables D, X and R are independent, 
so that (D^)„,fe, (X^)„,fe and {R!^)n,k are three independent double arrays of i.i.d. random 

variables. 

We of course realise that, in a convincing model, the random variables D, X and R should allow for 
some interaction (dependence), but assumption (A4) is made for simplicity. 
Let now 

k k 

Dn{k):=Y,Di and JJ„(A;) := ^ i?^„, (1) 

denote the total number of offspring and the total resources produced by generation n, respectively, 
given that generation n counts k individuals. The i.i.d. assumptions for random variables within the 
same double array allow us to use the shorter notations D{k) = Dn{k) and R{k) = i?„(fc) whenever 
wc limit our interest to their distributional prescriptions. Conversely, this is understood throughout 
the paper whenever we use this simplified notation. 

In general, the total resources produced by a generation may be insufficient to satisfy all the 
resource claims of the offspring. We define a policy as a function which determines a priority order 
among offspring, i.e. a rule to distribute among the next generation the resources produced by the 
current generation. 
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Definition 3.1 (Global definition of a policy). A policy is any function of the form 

TT : IJ ({t} X (]R+)*) ^ U {{t} X Sym(t)) : {t, (xfe)Li) ^ {t,Mi^k)i=i)), 

tGN teN 

where Sym(t) denotes the set of all permutations of [t] {1, . . . ,t}. 

In this definition, t corresponds to the number of offspring, and {xk)l.^i to their respective resource 
claims (i.e. independent realisations of the random variable X). The permutation '^t{{xk)k=i) € 
Sym(i) then gives the priority order that the society has chosen to satisfy the claims of the offspring: 
the individual T^t{{xk)l.^i){l) is the first served, etc. If s denotes the total of resources produced by 
the previous generation, the number of offspring having their claims completely satisfied thanks to 
the society's policy tt is thus defined by 

[max |1 <k<t: }_^^^^ ^Mixk)Ui)U) - ' otherwise. 
Note that this function necessarily satisfies 

Q-(O,0,s) =0-Q-(t,(xfc)Li,0) and < Q-(i, (xfc)Li, s) < 
for all s £ M+, aU t G N and aU {xk)l^-^ e (M+)*. 

We now recall (see Section 2) that all the offspring that are not completely satisfied, and only these, 
leave the society forever. This leads to the definition of the following stochastic process. 

Definition 3.2 (Global model). If tt is some policy, the resource dependent branching process 
(RDBP) on {D^,X^,Rl)^^k controlled by TT is defined as the integer-valued non-negative stochastic 
process (r„)„gN, defined by Fq = 1 and recursively 

r„+l = Q"" {Dn{Tn),(X^)^^^^"\Rn{Vn)), 

where Dn{-) and Rn{') are given by Equation (1). 

Such processes are Markov (see Proposition 5.1 hereafter). 



A remark on multiparameter policies 

According to our definition, a policy can only depend on the available resources and on the claims of 
the offspring. However, in more realistic models, the offspring could be characterised by many other 
diflFerent parameters, and it would be natural to allow a policy to depend on all these additional 
parameters. This is why, although we do not pursue such general models in this paper, we want to 
indicate shortly how we could adapt our definitions accordingly. 

We thus consider a new double array (y^)n,fc of random p-vectors defined on a corresponding 
probability space (17, J", P), such that all quadruplets (Z?^, X^, 1^, i?^) are i.i.d., with the same 
distribution as some quadruplet {D^X,Y,R). Here, the components of the random vectors 
(k G No) correspond to the different characteristic parameters of each individual (labeled by k) in 
generation n. 

Definition 3.3. A p-parameter policy is any function of the form 

TT : \J{{t} X (M+ X BY) ^ \Ji{t} X Syni(i)) : {t, ixk,yk)U,) ^ (t, 7rt((xfc, yfc)Li)), 

teN ten 

where p > 0, yk = {y'j}\ ■ ■ ■ , Uk'^), and B cM^ the set of possible parameter values. 
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The associated counting function is thus defined by 

[max |1 <k<t: Y.j=i ^M(x^,.m]U^){j) ^ ' otherwise. 

Definition 3.4. If tt is a p-parameter pohcy, the RDBP on (X^, y^, Z?^, i?^)fe_„ controlled by tt 
is defined as the integer-valued non-negative stochastic process (r„)„gN, defined by Fq = 1 and 
recursively 

r„+i = Q-(D„(r„),(x^itf)f:f"\i?„(r„)). 

Such processes are again Markov. As an example, we might think of some "Darwinian policy", 
which gives the priority to the claims of the (in several respects) fittest individuals. 



4 Particular policies 

In the following, we define policies we think of as being of particular interest. The first will be a 
neutral policy, which we call the first-come-first-served policy. It will serve as a point of comparison 
with the weakest-first policy and the strongest-first policy defined later. The latter two policies turn 
out to be extreme in a sense to be clarified. The reader will notice that, in a first approximation, 
these two extreme policies have some important features in common with the principles of certain 
society forms we believe to know rather well. 



4.1 First-come-first-served policy 

The fcfs-policy is a neutral policy in the vein that it serves the claims according to their respective 
arrival times. To exclude ambiguities in the definition, these arrival times of claims are seen as lying 
at the beginning of each generation, being almost surely different, and all preceding the times of 
producing offspring. 

Definition 4.1. The first-come-first-served policy (fcfs-policy) is the deterministic policy tt^ defined 
by^F((a;,)Li)=idM.^ 

The associated counting function C := Q'^ thus reads 



j 0, if i = or xi > s 
C{t,{xk)k=i,s) ^ < r fc I 

I sup <:l<K<t: i^j^i Xj < s> , otherwise. 

Definition 4.2. The first-come-first-served process (fcfs-process) on {X!^, D^^, R^)n.k is the RDBP 
controlled by tt^ , i.e. the stochastic process {Un)neN defined by Uq — 1, and recursively by 

Note that C{t, {X^)\^^, s) is a stopping time with respect to the natural filtration {J-e)i, where J-e 
denotes the cr-field generated by the X^'s for 1 < k < £. 

Interpretation and properties. The fcfs-society may be seen as a model of a laissez-faire society. 
Individuals are born, arrive at (a.s.) different times within their generation at maturity and submit 
their random resource claims. This continues as long as resources are available. Since the claims 
are i.i.d. random variables, it is not the society but the scarcity of resources which may impose 
constraints. This process has some similarity with the GWP, because, for given distributions of 
resource production and claims, the claims curtail the (effective) mean m of the intrinsic offspring 
distribution {pk)k- This is why, among the special processes we consider, the fcfs-process has the 
easiest structure. 

^The notation tt'^ should remind of the unordered x\,. . . ,x\ used in the definition. 
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4.2 Weakest-first poHcy 

The weakest-first policy (wf-pohcy) is an extreme pohcy, giving priority successively to the least 
demanding currently remaining offspring. 

Definition 4.3. The weakest-first policy (wf-policy) is the deterministic policy tt^ defined by 

'^V {{^k)k=i) = o'l where cr is the permutation of [t] such that x„(^i-f < . . . < .Tcr(t). 

Throughout this paper, for i.i.d. realisations {xkYj^^i random variable X, the increasing order 

w 

statistics will be denoted by xi^t < 2:2,* < • • • < xt^t- The associated counting function A'' := is 
now 

fO, if f = or .xi t > s, 



sup |l < A; < t : X)j=i ^j,t ^ i otherwise. 



Definition 4.4. The weakest-first process (wf-process) on [D^, X^, R!^)n^k is the RDBP controlled 
by TT^, i.e. the stochastic process {Wn)nGN defined by Wq = 1, and recursively by 

Note that A''(-, •, •) counts the maximal number of increasing order statistics of the random sam- 
ple {xk)l.^i which, starting with the smallest, can be summed up without exceeding s. Further, 
N{t, (X^)^^;^,s) is a stopping time on the filtration (J-"/)^ say, generated by the £ first increasing 
order statistics from all order statistics, beginning with the smallest one, but it is not a stopping 
time with respect to the natural filtration {Te)e. 

Interpretation and properties. The policy of the wf-society is to support always the weakest, fn that 
respect it comes close to the ideas of socialism and communism. In each generation, individuals are 
ordered according to their resource claims, and these order statistics are highly dependent of each 
other. Note that now a comparison with a GWP is (generation-wise) no longer possible. 

4.3 Strongest-first poHcy 

The strongest-first policy (sf-policy) is another extreme policy, giving successively priority to the 
most demanding currently remaining offspring. 

Definition 4.5. The strongest-first policy (sf-policy) is the deterministic policy defined by 
""f ((^A;)i;=i) = c, where a is the permutation of [t] such that 0:^(1) > • • • > x^^ty 

The associated counting function M := Q'^^ becomes 

f 0, if t = or .Tf f > s 
Mit,{xk)i=i,s)^ { j,^,^, ^ \ ,u ■ 

l^sup <.l<k<t: Y,j=t-k+i ^3,t < s> , otherwise, 

thus counting the maximal number of decreasing order statistics which can be summed up, starting 
with the biggest, without exceeding s. 

Definition 4.6. The strongest-first process (sf-process) on {X!^,D^,R'^)n,k is the RDBP controlled 
by n^, i.e. the stochastic process (<S'„)„£n defined by = 1, and recursively by 

= M(£>„(5„), {x!^)kd^-\Rn{Sn)). 

We note that M(t, (X^)^^-^, s) is a stopping time on the filtration generated by the first £ 

decreasing order statistics of all present claims, beginning with the largest one. It is again no stopping 
time on the natural filtration {Tt)t- 

Interpretation and properties. The sf-society is the model which serves the strongest individuals 
first. Recall that we identified the values of resource claims with the power to defend these claims. 
Hence this society shares important features with free-market policies and an uncontrolled capitalistic 
society. Since claims are again highly dependent, the technical difficulty in this model is comparable 
with the one evoked for the wf-society. However, as we shall see, there are non-negligible differences. 
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5 Main Results 



Before we present our results, we should, viewing motivation, draw attention to the fact that some 
of them are deeper than they might look. 

Here is one example. Let (r„)„ be a process living under an arbitrary policy tt. Since the sf- 
society is clearly the most restrictive one for the number of offspring which can stay, one feels that 
(r„)„ should always do at least as well as the process (S'„)„ governed by the sf-policy. Now, if for 
instance sf-process (<S'„)„ and (r„)„ have the same number k of individuals at time n, then it follows 
from the definition that r„_|_i is at least as large as Sn+i- Hence we expect in generation n + 2 on 
average more offspring from r„_|_i than from Sn+i- But then the extreme claims of the offspring 
of r„+i may be much larger than those from the offspring of Sn+i so that the inequality may, in 
generation n + 2, point to the opposite direction. 

Throughout this Section, all RDBPs are supposed to be controlled by some policy tt on some 
double array (D^, X^, Rn)n,k of i.i.d. triples of random variables satisfying assumptions (A1)-(A4) 
of Section 3. 

We begin with the easier results. 

It is important to first point out that any RDBP shares the following property which is so typical 
for most branching processes, namely either it explodes, or it becomes extinct. 

Proposition 5.1. Any RDBP (r„)„ is a Markov process. Moreover, it tends a.s. either to or to 
oo. The same result remains true in the multiparameter case. 

Proof. See Section 7.1. □ 

In accordance with Hypothesis 1 of Section 1, we must first answer the question under which 
conditions a given RDBP (Tn)n can survive, i.e. we must determine when the extinction probability 



lim r„ = 



Pn = 1 



is equal to 1, which will always be a central question in this paper. Note that, in the case when 
gr < Ij the probability of extinction can intuitively be made arbitrary small if we replace the initial 
setting Pq = 1 by Po = M for M sufficiently large. This seems hard to prove in all generality. For 
the wf-process, however, this is true and a trivial consequence of the following stronger result: 

Proposition 5.2. For all M > 0, Af G N, 



lim W„ = 



Wo = M 



Proof. See Section 7.1. □ 

A similar result can easily be shown for the fcfs-process. Nevertheless, determining a precise upper 
bound for P [lim„_>.oo Pn = | Pq = M] which goes to as M goes to infinity, remains an open 
question for the sf-process and general RDBPs. 



5.1 Multiparameter policies 

If the additional parameters of an individual are assumed to be independent of its number of offspring, 
its resource claim and its resource production, and if we consider some multiparameter policy that 
only depends on these additional parameters (but not on the resource claims), then the associated 
RDBP has exactly the same behaviour as the fcfs-process: 

Proposition 5.3. Assume that Y is independent of D, X and R, and that the multiparameter policy 
TT can be written under the form TTt{{xk,yk)k=i) — p{{yk)\=i)- Let {Tn)n denote the associated RDBP 
oniDtXtYj:,R>^)n,k. Then 

P[P„ m|Po = mo] ^ F[Un = m\Uo = mo] Vn,m,mo, 

where ([/„)„ is the corresponding fcfs-process. In particular, — Qu- 
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Proof. See Section 7.2. □ 

In general, the dependence may of course lead to highly complex situations, which we shall not 
study in this paper. Nevertheless, the wf-policy and the sf-policy being extreme policies, they give 
relevant sharp bounds for all general cases. 

5.2 General bounds 

It turns out that the wf-process is always an upper bound for any other RDBP, in a very strong 
sense. More precisely. 

Proposition 5.4. Let (r„)„ be any RDBP on {D^,X^,R^)n,k, and let (Wn)n be the wf-process on 
the same double array. Then, for all n, we have r„ < Wn a.s. In particular, qyy < q^. The same 
result remains true in the multiparameter case. 

Proof. See Section 7.3. □ 

We now turn to a comparison between (r„)„ and the corresponding sf-process {Sn)m which is a 
much more subtle problem. Indeed, we have already pointed out that it is in general not true that 
Sn < r„ a.s. for all n (see Section 7.3 for an explicit counterexample). Therefore all our attempts 
to compare trajectories turned out to be fruitless. 

We found it highly interesting that, nevertheless, we can prove that 

= 1 ^ gr = 1 ^ gs = 1 

(see Lemma 7.1 in Section 7.3), and remarkably, this suffices to deduce the following much stronger 
result: 

Theorem 5.1. Let (r„)„ be any RDBP on {D^, X!^, R'^)n,k, and let (W„)„ and {Sn)n be the wf- 
process and the sf-process on the same double array. Then, 

lim Sn < lim r„ < lim Wn a.s. 

In particular, qw l£ qr l£ Qs- The same results remain true in the multiparameter case. 

Proof. See Section 7.3. □ 

Such bounds are of considerable theoretical interest, and, as we shall now see, they are also serving 
as useful directives for individuals who have decided to adapt a specific type of society. Indeed, if the 
probability laws of the random variables {D^)n,k, {Xn)n,k and {Rn)n,k are fixed up to their mean 
m, /z and r, respectively, then it is in practice interesting to determine the critical mean resource 
production rr,c(™, A*): say, relative to the RDBP (r„)„, i.e. the value such that 

qr = I if r < rY.dm,^) and < 1 if t > rr,c{m,^j,). 

By Theorem 5.1, the following can be deduced: 

Corollary 5.1. For all m,fi, we have 

rwA^^l^) < ^r,c(?7T-,M) < Ts.ci'^^lA- 

Therefore, the study of the two extreme RDBPs gives non-trivial information about general RDBPs, 
without having to understand every single possible policy (see examples in Section 6). In the sequel, 
we shall thus essentially restrict our attention to the wf-policy and the sf-policy, as well as to the 
neutral fcfs-policy as a point of comparison. 

As we shall see in the next subsections, the computation of the critical mean resource production 
even shows more. The point is that the mean claim value plays only one part but that the resource 
claim distribution function F (which determines the mean, of course) plays itself an important part. 
Hence society may try to take influence on individuals to settle, under a fixed mean claim /i, for a 
distribution F which favours survival. 
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Remark 5.1. If D, X and R were not assumed to be independent, Theorem 5.1 would in general 
not remain true: the wf-policy and the sf-policy would a priori not remain extreme policies. We 
could then naturally wonder how different dependence patterns yield different extreme policies. Such 
questions may attract interest for further studies. 

5.3 The wf-process 

Theorem 5.2. Let {Wn)n be the wf-process on (_D^, X^, i?^)„_fc, where this double array of i.i.d. 
triples of random variables satisfies assumptions (A1)-(A4) of Section 3. Suppose m > \ and 
At > 0. 

(a) If r < mn and if t is the solution of 

rxdF{x) = -, (4) 
Jo m 



the 



(i) if mF{T) < 1, then qw = 1; 

(a) if mF{T) > 1 and F{r/k) > for some k >2 with p^ > 0, then qw < 1- 



(b) If r > mfj, and F{r/k) > for some k>2 with p^. > 0, then qyy < 1. 

Further, if there is no extinction, the process explodes a.s. and behaves asymptotically like a super- 
critical GWP with a new reproduction mean fh{> 1), say, defined by 



m — 



I, if r > mfi, 

iF(t), if r < nifj, and mF[T) > 1. 



Proof. See Section 7.4. □ 

The following remarks will provide a better understanding of these results. 
Remarks 5.1. 

i) We have assumed here that the resource claims are bounded (recall A3). However, we do not 
use this assumption in the proof of the above result. 

ii) The case m < 1 is trivial because then {Wn)n is stochastically smaller than a subcritical GWP 
and therefore bound to die out, since po > 0- Further, since all resource claims are non-negative, 
the case /i = is trivial too, because, if (X^)„^fc consists only of O's, then the process coincide 
with the standard GWP: in this case, survival is possible if and only if m > 1. 

iii) Among the non-trivial cases, the case (b) is the most intuitive one. Indeed, the condition 
r > m/i means that a typical ancestor produces in expectation more resources than his offspring 
will claim together. Consequently, when the population grows the law of large numbers ensures 
that the process will finally behave like a supercritical GWP, the asymptotic properties of which 
are well understood; see Bingham and Doncy [1974]. The additional condition po + Pi < 1 
implies for m > I that pk > for some fc > 2 so that the process can trigger off and reach 
any size with positive probability. This is a necessary condition for the preceding argument to 
hold; this condition becomes redundant if we replace the initial setting Wq = 1 by Wq — w for 
w sufficiently large. 

iv) Theorem -5.2 is sharp in the sense that mF{T) = 1 is the exact separation point between a.s. 
extinction and positive survival probability. However, unlike what occurs with GWPs, it is here 
not immediate to see under which conditions on the law {pk}k and F the critical case implies 
a.s. extinction. Note that, for fixed m and F, the parameter r = T{r/m) is increasing in r, so 
that the equation mF^r) = 1 defines a critical mean resource production rw,c ~ below which 
qw = 1 and above which qw < 1 provided that the additional condition involving F{r/k) > 
is satisfied. 
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We can further prove the following intuitive result. 

Corollary 5.2. Let {Wn)n be the wf-process on {D'^, X^, R'^Jn.kj where this double array of i.i.d. 
triples of random variables satisfies assumptions (A1)-(A4) of Section 3. Assume m > 1, fi > 0, 
and F{r/k) > for some k > 2 with > 0. Then, if /i < r, we have < 1. 

Proof. See Section 7.5. □ 
5.4 The sf-process 

We now present the extinction criterion for the sf-process. Since we deal here again with a process 
depending on the partial sum behaviour of order statistics - now on the sum of the largest ones - 
we expect analogies. Since is a lower bound for resource claims, the smallest order statistics are 
naturally bounded below. To facilitate a comparison between the sf-process and the wf-process we 
had made the assumption (recall A3) that resource claims are bounded above, i.e. that F has a 
finite support. This assumption facilitates the mathematical treatment, but is also highly reasonable 
as far as the interpretation of the model. 

Theorem 5.3. Let {Sn)n be the sf-process on {D'^, X^, R^)n,k) where this double array of i.i.d. 
triples of random variables satisfies assumptions (A1)-(A4) of Section 3. Suppose m > \ and 
/X > 0. 

(a) Lf r < mjjL and if 9 is the solution of 

rb ^ 

xdF(x) = —, (5) 
m 



then 



(i) z/m(l - F{e)) < I, then qs = 1; 

(ii) if m{l — F{6)) > 1 and F{r/k) > for some k>2 with pk > 0, then qs < 1. 



(b) Lf r > mfj, and F{r/k) > for some k>2 with pk > 0, then qs <1. 

Further, if there is no extinction, the process explodes a.s. and behaves asymptotically like a super- 
critical GWP with a new reproduction mean m{> I), say, defined by 

m, if r > mil, 

m{l-F{9)), tfr<miiandm{l-F{9))>l. 
Proof. See Section 7.4. □ 
Remarks 5.2. 

i) As noted for wf-processs, the case to < 1 is trivial, implying qs = 1. The case /i = is trivial 
too: the process then coincides with the standard GWP, so that survival is possible if and only 
if TO > 1. 

ii) Again, the additional condition involving F{r/k) > just serves to ensure that the process 
can grow. It becomes superfluous if we replace the initial setting 5*0 = 1 by = s for some s 
sufficiently large. 

iii) The critical case is now determined by the equation to,(1 — F{9)) — 1. Note that, for fixed 
TO, the parameter 9 — 9(r/m) is decreasing in r, in the same way that, in Theorem (5.2), 
T{r/m) was increasing in r. Again, the equation to(1 — F{9)) = 1 thus defines the critical 
mean resource production rs^c- 

We have, and here intuition is correct, the following result. 

Corollary 5.3. Let (S'„)„ be the sf-process on (_D^, A^, _R^)„ fc, where, again, this double array of 
i.i.d. triples of random variables satisfies assumptions (A1)-(A4) of Section 3. Assume to > 1, 
/X > 0, and F having finite support [a,b], < a < b < oo. Then, if r < fi, we have qs = 1. 

Proof. See Section 7.5. □ 
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5.5 The fcfs-process 

As a term of comparison, it is interesting to observe what happens in the case of a fcfs-process. 

Proposition 5.5. Let (J7„)n be the fcfs-process on {D^, X^^, Rj^)n^k, where this double array of i.i.d. 
triples of random variables satisfies assumptions (A1)-(A4) of Section 3. Suppose m > 1 and ^ > 0. 

(a) Ifr<fi, then qjj — 1. 

(b) Ifr>fj,, then qjj < 1. 

Further, if there is no extinction, the process explodes a.s. and behaves asymptotically like a super- 
critical GWP with reproduction mean m. 

Proof. This foUows immediately from the well-known properties of the GWP and the strong law of 
large numbers. □ 

Remarks 5.3. 

i) We have assumed here that the resource claims are bounded (recall A3). However, we do not 
use this assumption in the proof of the above result. 

ii) As noted for wf-processes, the case m < 1 is trivial, implying qu — 1. The case /i = is trivial 
too: the process then coincides with the standard GWP, so that survival is possible if and only 
if m > 1. 

iii) The critical mean resource production is now simply defined by rjj^c — 



6 Examples 

We now give examples. It will be interesting to notice that the critical mean resource production for 
a wf-process turns out to be lower than one would intuitively expect, and conversely for a sf-process. 

i) Let F be the uniform distribution function on {0,d), say. Then /j, = d/2. As in Theorems 5.2(a) 
and 5.3(a), let r < nifi = md/2 and suppose that F{r/k) > for some k > 2 with > 0. 
First, focus on the corresponding wf-process. The value of r (see Equation (4)) is thus deter- 
mined by: 

r 



[ xdF{x) — [ -xdx 
Jo Jo d 



2d' 



which yields t = ^2dr jm. Therefore, Fij) — \/2r j (md). The critical mean resource produc- 
tion r^o c is thus determined by 



mF(r) = 1 ^ y2r^.cm/d=l 

which implies r^^c = d/2m = fx/m. Note that, for larger m, r^.c is quite close to the smallest or 
second smallest (expected) claim of offspring of two families. Indeed, r^^c — d/2m > d/ {2m-\-l), 
which is the expected value of the smallest order statistic of 2m i.i.d. U[0, d] random variables, 
and r^j^c < 2d/{2m -\- 1). 

With such a low productivity of resources, the sf-process and the fcfs-process would die out 
very quickly, as we shall see now. 

In the case of a sf-process, we need to determine 9, which is defined by (see Equation (5)): 



r 



1 , d^-e^ 

-xdx = — — - — 
mod 2d 



and so, straightforward calculations yield rg^c — d{l — l/2m) = /i(2 — 1/m). 

We note that the critical mean resource production is now 2m — 1 times higher than for the 
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corresponding wf-process. Moreover, comparing with the corresponding critical mean resource 
production for a fcfs-process, r„ ^ = M; gives 

ru,c - rni,c = m(1 - = r^^c - Tu.c- 

Figure 1 compare the behaviour of r^_c and r^^c as functions of m. The area between the two 
curves corresponds to a critical area, where the population can both survive or die, depending 
on the policy. This shows how the study of the two extreme RDBPs yields highly non-trivial 
information about general RDBPs, without having to understand every single possible policy, 
as already pointed out in Section 5.2. 



Figure 1: For /i = 1, the critical mean resource productions r^^c and r^ c are plotted (in blue and 
in red, respectively), as functions of m. 



ii) Of course, we realise that the uniform distribution pushes the largest and the smallest order 
statistics far apart, namely, for n observations, by a factor n (in expectation). Therefore, we 
find it informative to look at the case when the resource claim distribution is more concentrated 
around its mean, as for instance in the case of a beta distribution (on (0, 1), with parameters 
a and b). 

The distribution function is then defined on (0,1) by the regularised incomplete beta function: 
F{x) = Ia,b{x). The mean resource claim is given by /i = As in Theorems 5.2(a) and 

5.3(a), let r < m/i = and suppose that F{r/k) > for some k > 2 with pk > 0. 

First, focus on the corresponding wf-process. The value of r (see Equation (4)) is determined 
by: 

r ^Pf ^ r ^"(l-^)"'^ B{a+l,b) ^ a 
m Jq Jq B{a,o) B(a,b) a + o 

which yields r = la+i b { — ^-^). The critical mean resource production r^, c is thus determined 
by 

a + b'" 



mFir) = 1 ^ mla.b [Jalub (^^-.-^ j j " ^' 
which implies 

QjTTl , -I , , , , 

ru,,c = ^^^4+i,fc(4,h(l/"i))- 
Now look at the corresponding sf-process. The value of Q (see Equation (5)) is determined by: 

-=/ xdF(a;) = ^(l-/,+i,b(f?)), 

which yields B ~ ^a+i b{^ ^ m^^^)- Straightforward calculations then give the critical mean 
resource production r^^c = ^(1 - la+iA^aM^ ^ V"i)))- 
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Observing that Iaj{l — x) — 1 — /^^^(a;), we deduce that }3(1 — a;) = 1 — ^(x). The formula 
for c can thus be rewritten as 

am / -1 / \\ 

For a fcfs-process, the corresponding critical mean resource production simply reads r„ ^ = A*- 
Further, for big to, we could use the approximation 

i.A^) = (- + + oiz')) 

B{a, p) \a a + 1 ) 
(see, e.g., Pearson [1968]). Straightforward calculations then give, at leading order: 

r^.c - (-Bia, b)) + 0(m-2/-), 

a + 1 Vto / 

and 

b / b \^^^ 

Ts.c = 1 + 7— (& - a{b + 1)) -B{b, a) + ©(to-^/"). 

0+1 \TO / 

Figure 2 shows the critical areas in some typical cases, as the peak is centered, moved to the 
left or to the right. Figure 3 shows, in the centered case, how the critical area narrows as the 
dispersion around the peak diminishes. 




Figure 2: The critical mean resource productions r^^c and r^^c are plotted as functions of to, in the 
case of a B{a,b) resource claim distribution, for typical values of (a, 6): (2,10) in blue, (14,14) in 
yellow, and (10, 2) in red. 



iii) Let F be the distribution function of an exponential random variable with parameter A. This 
is not bounded and our results in the sf-case cannot be used. However, it could be interesting 
to see what happens for the corresponding wf-process. The mean resource claim is given by 
fj, = 1/X. As in Theorem 5.2(a), let r < to/i — m/X and suppose that F{r/k) > for some 
k>2 with pk > 0. 

The value of t is determined by (see Equation (4)) 

- = / xdF(x) = / Xxe-^^'dx ^l-efr+l-] g-^f^+i/^) 
m Jo Jo X \ XJ 

which yields t = —j (l + W [— — ^)]), where W[-] denotes the Lambert W function (see, 
e.g., Corless, Gonnet, Hare, Jeffrey, and Knuth [1996]). The critical mean resource production. 
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Figure 3: The critical mean resource productions r^j^c and r^^c are plotted as functions of m, in the 
case of a B{a,b) resource claim distribution, for different symmetric values of (a, &): (1, 1) in blue, 
(2, 2) in pink, (3, 3) in red, (4, 4) in orange and (10, 10) in red. 



r^i.c, is thus determined by 

mF{T) = 1 m — em exp 

After simplifications, we get r^^c = y (^1 ^ ('rt — 1) log (^^^f^jj j , where we recall that ^ — fi. 
For big TO, this becomes r^o ^ w 1/Ato = fi/m. 



X I 

e \X 



' w,c 

m 



= 1. 



7 Proofs 

7.1 Preliminary results 

We prove here Propositions 5.1 and 5.2. 

Proof of Proposition 5.1. Assume tt to be some policy. Given r„, the distributions of Z?„(r„) and 
i?„(r„) are independent from ri,...,r„_i, so that r„+i = Q''(Ai(r„), i?„(r„)) is 

independent of Fi, . . . ,F„_i, given F„. Thus, (F„)„ is a Markov process. 

Now note that, since Q'"(0, 0, s) = and L'„(0) = for all n e N, we have {F„ = 0} C {r„+i = 0} 
so that is an absorbing state for the process (F„)„. Moreover, since 

F„+i = g'^p„(F„), {X'^)kZT\Rn{Tn)) < A.(F„), (6) 

it follows that 

P[F„+i = 0|r„] > P[i?„(F„) = 0|F„] = pI-, (7) 

where the last equality holds because of the assumption of independent reproduction. Therefore, the 
absorbing state is accessible from any state s € N, with at least probability p^ > 0. The state is 
thus the only absorbing state, and, as (F„)„ is a Markov process, we conclude 

P[0 < F„ < s i.o.] =0, Vs e No. (8) 

The same arguments immediately adapt to multiparameter policies. □ 
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We now turn to the proof of Proposition 5.2. 

Proof of Proposition 5.2. Let M e Nq and let Wn\ ■ ■ ■ , Wn^'^^ be M i.i.d. copies of a weakest-first 
process. We then have the following (superadditivity-type) inequality, namely, for all n, fc G No, 

P[W„ <k\Wo = M] < P[Ty,W + . . . + W^f ^) < k I W^^^ = ... = W^^^^ = 1], (9) 

which we shall prove first. 

To see this, we begin with the case n — 1. The Ihs of (9) becomes, by an additional conditioning 
on D^{M) = D{M), 

P[Wi < k\Wo = M] = P[D{M) <k]+ P[Wi < k\D{M) > k,Wo ^ M]P[D{M) > k], 

since the distribution of D{M) depends only on M, and the probability that Wi does not exceed 
D{Wo) = D{M) equals 1. 

Now suppose we do the same conditioning on the rhs of (9), that is, for the offspring of the M 
partitioned processes. The first term P[D{M) < k] is the same on both sides, as reproduction is 
independent. Hence we can limit our interest to the corresponding second term with more than k 
offspring. 

The distributions of the total created resource space and of the claims are also by definition the 
same on both sides; therefore it suffices to look for the moment at the influence of the order statistics 
of claims in a fixed sequence of claims on a fixed resource space R, say. In the Ihs-model the resource 
space R is global (i.e., united) whereas in the rhs-model it is local (i.e., compartmented). On the 
Ihs the count of individuals to stay is the count of the globally smallest order statistics of claims 
which can be accommodated by R whereas on the rhs the count is on the locally smallest order 
statistics of claims. The latter, put in increasing order, are a subsequence of the sequence of claims 
in increasing order. Hence the rhs count cannot exceed the Ihs count. Passing from the count to the 
corresponding probability measures on both sides proves (9) for n = 1, that is, Wi is stochastically 
larger than W^^^ + ... + W^'^^\ 

But now, the inequality (9) must hold in particular if we replace on the rhs the number M by 
some M' with M' < M. Hence the stochastic order is maintained through the next generation, and 
thus, by recurrence, through all generations. This implies that (9) is true for all n S No. 

Finally, choosing fc = in (9) and taking the limit for n — >■ oo, we obtain by independence of the 
processes {Wn^)n that 

P[W„ ^ I M^o = M] < P[wi^^ ^ 0, . . . , Wj,"'^ ^ I W^o'^' = . . . = W^"-'^ - 1] = q^, (10) 
which completes the proof. □ 

Remark: The superadditivity-type inequality (9) is no more correct if the wf-process is replaced 
by the sf-process. Indeed, a very large claim may now force on the rhs all the offspring of one 
subpopulation to leave, but this effect stays still local whereas it may be large on the global Ihs. 
The same argument implies that the corresponding inequality of (9) is no longer true even for 
the fcfs-process. This exemplifies at the same time the adherent difficulty in estimating extinction 
probabilities for arbitrary policies. 

7.2 Multiparameter policies 

We now prove Proposition 5.3. 

Proof of Proposition 5.3. We prove by induction on n that 

P[r„ = m|ro = mo] = P[C/„ = m|[/o = mo] 
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for all n, m, tuq. It is clear for n = 0. Moreover, if it is true for some n, then we can write 
P[r„+i = r7i|ro = Too] 

nOO POO 

= E E / ■ ■ • / d^(^i) ■ ■ ■ di^(^t)P[pt((n)Li) = f^lPfr™ = fc|ro = TOo]P[D„(fc) = <] 

P[r„+i = m|ro = mo,r„ = k,Dn{k) ^ t, ^ xi, . . . , ^ xt, pt((Yfe)^^i) = cr] 
= 5Z E / ■■■/ <iF{xi)...dF{xt)^P[Un = k\Uo = m„]P[D,,{k)=t] 
?[[/„+! = m\Uo = Too, C/„ = k,Dn{k) = t, X,^, = . . . ,X* = Xa.(t)] 

OO 

= J2 = ^I'^O = TOo]P[A.(fc) - t]P[Un+i - m|C/o - TOO, Un = fc, D„{k) = i] 

fe,t=0 

=f'[C^n+i = m\Uo = Too], 

and thus the statement is proved. □ 

7.3 General bounds 

We first prove Proposition 5.4. 

Proof of Proposition 5.4- Let tt be any policy (the same arguments immediately adapt to multipa- 
rameter policies). First note that, by definition of N{-) and M{-), 

M{t, (xfe)Li, < QHt, {xu)U,,s) < N{t, {xk)U„s), Vi, V(xfc)Li, Vs. (11) 

We shall now show by induction that if (r„)„ is the RDBP controlled by tt, it follows that r„ < Wn 
a.s. for all n, given Wq = Fo = ^o = 1. Indeed, it is true at any time at which the three processes 
have the same number of individuals, and hence for n — 0. Now, if it is true for some n, we deduce 
that, a.s., 

r„+i - Q"(A.(r„), (x^)t"r"\i?„(r„)) (12) 

< iv(A.(r„), (x„^)t"r"\i?„(F„)) (13) 

< N{DjWn), {Xi:)^2r-\Rn{Wn)) = (14) 

as the mapping (i, s) N{t, {xk)l.^i, s) is increasing in both arguments. Hence the inequality is 
also true for n + 1. 

It follows that P(F„ < Wn) = 1 for all n. Since the limiting extinction probabilities of (F„)„ and 
(W^n)n must exist, we must also have qw < gp- D 

We now give an explicit counterexample showing that it is in general not true that Sn < F„ a.s. 
for all n, given 5o = Fg = 1. The underlying idea was already explained at the very beginning of 
Section 5. 

Counterexample. We have assumed pk > for some k > 2 (see (Al)); to fix ideas, assume 
that p3 > (the argument can be adapted in any case). Then consider the deterministic policy tt 
given by 

V(3), if j = 1 and i > 3, 

,f I cr(l), if 7 = 2 and i > 3, 

MMUiKj) ^ { „ ' ./, „ , - ' (15) 
(t(2), if j = 3 and t > 3, 

^cr(j), otherwise. 
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where a is the permutation such that x^^^) ^ ■ ■ ■ ^ Xa(t) (i-e-, by definition, a = Ti'f ((a;^)^^]^)). Let 
for example 

Dq = 3,Xq > Xq > Xq,Xq + Xq < i?Q < + Xq 



and then 



Dl =Dl = i, Xl + Xl + Xl < R\ 

v& ^ Tjl I d2 



These events will occur simultaneously with positive probability, as ps > 0. But then we immediately 
see that r2 = < 3 = 5*2 in this case. 

We now turn to the proof of Theorem 5.1. For that purpose, we shall need the following interesting 
result, which we prove first. 

Theorem 7.1. Let (r„)„ be any RDBP controlled by some policy tt on {D'^, X^, R!^)n.k, where this 
double array of i.i.d. triples of random variables satisfies assumptions (A1)-(A4) of Section 3. Let 
1 > > 0. For all 5 > 0, there exists some Lg £ No such that 



Tri+l ^ Sn+1 r 



r„ > Lf., Sn > La 



> 1 - ??• 



The same result remains true in the multiparameter case. 

Proof. Let 1 > > 0. Since j:D{k) — ^ m a.s. and \R{k) — ?> r a.s., we have, for all e > 0, 









sup 


^D(k) - m 
k 


< e 


_k'>m 





1 and P sup rR{k) — r < e 1, 

_k>m k 

as m — >■ oo. Therefore, for all e > 0, we can find E No such that 

{m - e)k < D{k) < (m + e)k and {r - e)k < R{k) < {r + e)k 



(16) 



(17) 



must hold (simultaneously) for all k > K^, with at least probability 1 — ry/2. 

On the one hand, if we assume that J2k=i -^n ^ Rik) is given, we necessarily have, given Sn ~ k, 

M(L(m - e)k\ , (X,t)Er^^'^ , Rik)) < Sn+i < A/(r(m + e)k^ , (X^)[.LT'^'^ > m), (18) 

where [x\ denotes the largest integer < x and where \x~\ denotes the smallest integer > x. Since 
M{t, {X!^)j,^^, s) is stochastically increasing in s, we thus have 

M(L(m - e)k\ , (X,t)iLT'^'^ > (r ~ e)k) < Sn+i < Mi\{m + e)k^ , (X,t)iLT'^'^ - + ^)k) (19) 

for all k > K^, with at least probability 1 — 7//2. Now, using this and Theorem 2.3 of Bruss and 
Robertson [1991], we deduce that, for all S > 0, there exists some Kg ^ > such that, for all 
k> K', 



5,ei 



[(TO-e)fcJ [l-F 



(r — e)k 



[{m-e)k\ 
< Sn+1 < \{m + e)k'] (l-F 



- (5/ (2m) 
(r + e)k 



\{m + e)k] ' ' 

with at least probability (1 — ?7/2)^ > 1 — where the continuous function 9{-) is defined by 

rb 



xdF{x) 



(20) 
(21) 

(22) 
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Taking e = > small enough yields on the one hand 

mfc (1 - F {e{r/m)) ~ 5/m) < 5„+i < mfc(l - F{9{r/m)) + S/m) (23) 

with at least probability 1 — ?7, for all k > Kg =: Kg. 

On the other hand, if we assume that J2f=i < R{k) is given, then, conditioned on the event 
{Sn — k}, we must have 

Mi\im + e)k^,{X^^i^:l+^^'\R{k)) < 5„+i < M(L(m + 6)fcJ,(X„^)L(™+'^)'=J,i?(fc)) (24) 

with at least probability 1 — 77/2, for all k > K^. 

A similar argument as above then shows that, for all S > 0, there exists some Kg E Nq such that 

mfc(l - F{e{r/m)) - 5/m) < Sn+i < mk{l - F{e{r/m)) + 5/m), (25) 
with at least probability I — rj, for all k > Kg. 

Now let Ls :— \j\ax.{K'g, Kg). We have thus proven that, for all 5 > 0, given Sn — k, 

mk{l - F{e{r/m)) - 5/m) < 5„+i < mfc(l - F{9{r/m)) + 5/m) (26) 

must hold with at least probability 1 — rj, for all k > Ls. 
Note that, given — k > Lg, the inequalities 

r„+i = Q^Dn{k),{xX=['\Rnik)) (27) 
>M(i^„(fc),(X^)tf \i?„(fc)) (28) 
>mk{l-F{e{r/m))-5/m) (29) 

hold with at least probability l—rj. Combining Equations (26) and (29) then implies that, conditioned 
on r„ > Ls and Sn > Ls, the inequalites 

> m(l - F{9{r/m)) - 5/m) > - 25 (30) 

must hold with at least probability l — rj. This proves the statement of Theorem 7.1. □ 

We can now deduce the following lemma, which will be the key for Theorem 5.1. 

Lemma 7.1. Let (r„)„ be any RDBF controlled by some policy tt on (Dfj, X^', i?^)„jt, where this 
double array of i.i.d. triples of random variables satisfies assumptions (A1)-(A4) of Section 3. 
Then, if qr = 1, we must have qs ^ 1 too. The same result remains true in the multivariate case. 

Proof. By Theorem 5.3, we know that, given Sn — > 00, {Sn)n behaves asymptotically like a super- 
critical GWP with some reproduction mean m > 1. Therefore, Sn+i/Sn — >■ rh a.s. Let 1 > 77 > 0. 
We deduce that, for all (5 > 0, there exists some L'g e No such that, given 5*^ > L'g, Sk+i/Sk > m — 5 
must hold with at least probability 1 — 77. Combining this with Theorem 7.1, and putting 5 := m/A, 
we get, given S'n — >■ cx) and 5*^ > max(L5,L^) =: L, 

^^>^-5>m^25 = m/2 (31) 
1 fc Sk 

with at least probability (1 — 77)^ > 1 — 277. 
Define, for fc e N, 

._ finf{A^ : Sn>k\/n> N), if 3N such that Sn >k\/n>N, ^^^^ 
1+00, otherwise. 
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We can thus write 

P[r„ ^ cx)|ro - 1] > P[^„ ^ oo\So - l]P[r„ oo|S„ -> oo, 5o = To = 1] (33) 

oo 

- Y,P[^n^^\Sn^^,So^To^l,TL^N]P[TL=N\Sr,^^,So = l] (34) 

•P[5„^oo|S'o = l] 

> 51 E ^ ^IS-n ^,5o = Po - l,ri = N,Tm > L] (35) 

N M>N 

■ P[rM > L\Sn ^ OO, 5o = Po = l]P[rL = iV|5„ ^ 5o = 1]P[5„ ^ c»|5o = 1] 

> (1 - 2?7)P[^„ ^ oo|5o = MJ2J2 PPA^ ^ ^1^" ^ c^, 5*0 = Po = 1] (36) 

N M>N 

•P[ri=iV|^„^oo,5o = l] 
= (1 - 27y)P[^„ ^ oo|5o = 1] 5] PpA/ > -> ^, S-o = Po - 1] (37) 

M 

■ P[Sn >L\Jn> M\Sn -> oo, 5o = 1]. 
Now note that P[S'„ > L Vn > MIS',! -> oo, 5o = 1] ^ 1 as M -> oo. Let A/q G Nq such that 

P[5„ > L Vn > MQ\Sn ^ oo, 5*0 = 1] > 1 - r?. (38) 
We deduce, as this sequence is increasing in M, 

P[Pm > L\Sn ^ oo, 5*0 = Po = 1]P[S'„ > i Vn > M|5„ ^00,50 = 1] (39) 

> (1-7?) J2 P[rM>i|^„^(X),5o-Po = l] (40) 

M>AIo 

> (1 - 7?)P[3Af > Afo : Pm > i|^„ ^ 00, = Tq = 1], (41) 
so that finally, 

P[P„ oo|Po = 1] (42) 
> (1 - 3ry)P[3A/ > Mo : P^ > L\S„ ^ ^, 5o = Po = 1]P[5„ ^ ^|5o = 1]. (43) 

Given 5„ 00, it is clear that the process (Pn)n niust be able to reach or exceed any finite level 
with positive probability, so that 

P[3M > Mo : Tm > L\Sn ^ 00, So = Pq = 1] > 0. 

This and Inequality (43) allow us to conclude that, if P[Pn — )• oo|Po = 1] = 0, then P[5„ oo\So — 
1] = too. The conclusion then follows from our Proposition 5.1. □ 

We can now give a proof of Theorem 5.1. 

Proof of Theorem 5.1. To simplify notations, write W := lim„ Wn, P '.— lim„ P„ and S := lim„ S'„. 
As these three random variables take a.s. only the values and 00 by Proposition 5.1, we get from 
conditioning on P, 

P[S^ < P < W] = qrP[S < P < W\T = 0] + (1 - qr)P[S < P < W\T = 00] (44) 
= qrP[S = 0|r = 0] + (1 - qr)P[W = oo|r = 00] (45) 
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Now, by Lemma 7.1, qr — 1 implies 95 = 1, and thus F a.s. implies S — a.s. Further, by 
Proposition 5.4, we must have W — 00 a.s. if F = cx) a.s. Consequently, the preceding probability 
simply becomes 

P[S <T <W]=qr + {l-qr)^l, (46) 
thus proving the result. □ 



7.4 Proof of the extinction criteria 

In this section, we give a proof of Theorems 5.2 and 5.3. We shall make repeatedly use of the 
following lemma, which we shall prove first. 

Lemma 7.2. Let Xi,X2, . . . be i.i.d. real-valued non-negative random variables with mean /z < 00 
and continuous distribution function F. Further, let (<&n)n be a sequence of integer-valued random 
variables with $„ — > 00 a.s. as n ^ 00, and let {'^n)n be a sequence of real random variables with 
— >■ 00 a.s. as n — > 00. Suppose that ^'„/$„ — > p a.s. with < p < p, and that r is the solution 

of 

/ xdF{x) = p. 
Jo 

Let N{; ; •) be defined by (2). Then Ar($„, {Xk)tZi, ^ as- 



Proof Since ^„/<l>„ p > a.s. as n — >■ 00, we have 



Ve > 



sup -p| < e 



n>m 



1 a.s. as m — )■ 00. (47) 



To simplify notations, we write N{t, s) :— N{t, (-'i^fc)fc=i, s)- As this function is stochastically increas- 
ing in s, we deduce that, for all < e < p and all (5 > 0, the inequalities 

JV($„, (p - 6)<1>„) ^ ^($„,VI'„) ^ jV(j>„,(p + €)$„) 

$ - $ - $ ^ ^ 

must hold (simultaneously), for all n sufficiently large, with probability at least 1 — S. Refining a 
result of Coffman, Flatto, and Weber [1987], Bruss and Robertson [1991] (see Theorem 2.2, p. 615) 
have shown that 

^^lilhlA ^ F{t{s)) a.s. as n ^ cx) (49) 



where t(s) is the solution of 



/ xdF{x) = lim — =: s, (50) 

Jq n-i-oo n 



provided that the latter limit exists and satisfies < s < p = E[X]. Note that t(-) is of course 
continuous on (0,/i). 

As (f>„ — 00 a.s., the left-hand side variable of (48) must converge a.s. to i^(r(p— e)) and the right- 
hand side variable a.s. to F{t{p -f e)). Since e > is arbitrary and lim£_j.o+ F{t{p ± e)) — F{r{p)) 
by continuity of F{-) and t(-), the Lemma is proved. □ 



We can now prove Theorem 5.2. 
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Proof of Theorem 5.2. We first prove statement (a). Suppose r < mfi and mF{T) < 1, and let 

W'^ifl) = {u}efl: Wniuj) cx) as n ^ cx)}. (51) 

In the following, we write iV„(t, s) := N{t, {X^)l^j^,s). 
Now look at 

E[W„+i|W^„ =w]= E[Nr.{DJw),Rn{w))\Wn = w]. (52) 

Since Rn{w)/w — )■ r a.s. as — oo and I?„(w)/w — )■ m a.s. as w — > cx) and m > 0, we have 
Rn{w) / Dn{w) — 7- p = r/m a.s. According to Lemma 4.1 of Bruss and Robertson [1991], there exists 
a sequence T^j ^ t a.s. as ■«;—>■ oo with 

E[7V„(i?„(w),i?„(7«))|Ty„ ^w]< E[Dn{w)]FiT^) (53) 

where r is the solution of JJ" x dF{x) = p. Since F is continuous we can find, for each e > 0, a value 
w — w{e) such that F{Ty) < F(t + e) for all v > w. Thus, from equations (52) and (53): 

E[Wn+i\Wn^v]<mvF{T + e), v>w (54) 

and so 

oo 

E[Wn+l\Wn >W]^Y1 P[Wn = v\Wn > w] E[Wn+l\Wn = v] (55) 

V — W 

oo 

< mF(r + e) ^ wP[T4^„ = t;|VK„ > w] (56) 

V—W 

= mF{T + e)E[W„\Wn>w]. (57) 

Since mF{T) < 1, we can choose, again by continuity of i^, e > such that mF(r + e) < 1. The latter 
implies then that E[W„] must be bounded. Consequently, Proposition 5.1 allows us to conclude that 
P[M^°°(ri)] = 0, or equivalently qw — 1- This proves the first part of Theorem 5.2(a). 

To see the second part of Theorem 5.2(a), we now suppose that mF(r) > 1 and that F{r/k) > 
for some k > 2 with pk > 0. First note that under these conditions each individual can (alone) 
provide, with positive probability, enough resources for at least two offspring to stay. Therefore 
P[Wn > 2""-'^] > 0, so that {Wn)n can reach or exceed any finite state w with positive probability. 
Therefore it suffices to show that 

3a > 1 : liminf P[Wn+k > a''w\Wn >w]>0. (58) 

/c— >-oo 

Let h{j,a,w) := P[W^„4.j < a^w\Wn+j-i > a^^^w]. Then clearly 

P[Wn+k > a''w\Wn >w]>il- h{k, a, w))P[Wn+k-i > a''-^w\Wn > w] (59) 
and so by recurrence 

k 

P[Wn+k > a''w\W,, >w]> Y[{1 - h{j,a,w)). (60) 

i=i 

Therefore a sufficient condition for (58) to hold is 

oo 

h{j, g, if;) < oo and h{j,a,w) < 1 for all j. (61) 
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The second condition of (Gl) is clearly satisfied. We now show that the first condition also holds. 
Since Nn{-, ■), Dn+j{-) and R„^j(-) are stochastically increasing in their arguments, we have 

hij, a, w) = P[Wn+j < a^w\Wn+j^i > a^-'w] (62) 
= P[iY„+, (W„+j_i),i?„+,(T^„+j_i)) < a^w\Wn+j-i > a^'^w] (63) 
< F[N„+,{Dn+,{la^-'w\),Rn+j{la'-'w\)) < a^w]. (64) 

Choose e > such that (m — e)F{T) > 1 and put a — {m — e)F{T). Then, < 1 so that 



3 = 1 



= ;^P[D„+,(K-'u.J) = 0] <oo. 



(65) 



Since the Dn+j are independent random variables, it follows from the Borel-Cantelli Lemma that 
I'[Dn+j{[a^~^ w\) — i.o.] = 0. Therefore, for j sufficiently large, inequality (64) is equivalent to 



where 



and 



h{j,a,w) < P[Lj < Rj] 

Nn+j {Dn+o ( K^^H ) , Rn+j ( [a^~^w\ )) 



(66) 



(67) 



R,= 



Put Pj — [a-'^^wj and look first at the random variables Dn+j{Pj) / ^ Recall that 

D{n) is a sum of n i.i.d. random variables with finite mean m and finite variance af^, say. Also, 
since a > 1, /3j/j — > (X) as j — > oo. Therefore, it follows from the Hsu-Robbins theorem of complete 
convergence (see Theorem 1 of Hsu and Robbins [1947], and for furthergoing results Asmusscn and 
Kurtz [1980]) that D[[ij)/ j3j -> m completely, and thus 



V(5 > 



D{(i,] 



>S 



< oo. 



(69) 



Further, since a^w//3j — )■ a as j — > oo and _D„+j(-) ^ D{-), we obtain from (68) and (69) 



V(5 > : 



OG 



Rj 



a 
m 



> 5 



< oo. 



(70) 



Secondly, to study the convergence of Lj defined in (67) we turn to Lemma 7.2 with = D{/3j) and 
= R{f3j). Since ^j/Pj — m completely and '^j/Pj — r completely (again by the Hsu-Robbins 
theorem), we have 



p = — completely, as j — > oo. 



(71) 



Therefore, in particular, ^'j/'f'j p a.s., so that the conditions of Lemma 7.2 are satisfied. It follows 
that if Lj in (67) allows for a limit (in some sense) say, then we must have I = F(t), where r 
is defined as in Lemma 7.2. Using this and the Chernoff-type estimates obtained by Coffman ct al. 
[1987] with a = jS, we obtain after some straightforward simplifications: 



7V(,.) 



-F{r) 



> S 



< 2e ^'^M , 



(72) 
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Again A^^+j(-, ■) ~ A^(-, •) and Pj/ j — )■ oo, and thus Lj — )■ F{t) completely as j — )■ cx:, implying 

oo 

V(5>0: ^V[\Lj-F{t)\> 5]<oo. (73) 

Now choose 5 = ^\F{t) — a/m\ > 0. Note that the event {Lj < Rj} cannot happen unless 
\Lj - F{t)\ > (5 or \Rj - a/m\ > 5. Therefore, from (66), 

hU,a,w) <P[L, <Rj] (74) 
< P[|ij - i^(r)| > (5] +P Rj-^ >S (75) 
so that, accordmg to (70) and (73) 

oo 

/i(j, a, w) < cx). (76) 

This completes the proof of statement (a). Statement (b) is easily obtained, using Theorem 2.1 
of Bruss and Robertson [1991]. □ 

Theorem 5.3 can be obtained by similar considerations, the role of r being now played by 0, 
defined by Jg xdF{x) = ^, as in Theorem 2.3 of Bruss and Robertson [1991]. 

7.5 Proof of Corollaries 5.2 and 5.3 
Lemma 7.3. Assume r > ii and let r be defined by 

xdFix) = — . 

m 

Then, mF{T) > 1. 

Proof. As F{x) > for all x > a and as t > a (because xdF{x) — ^ > > 0), we deduce 



, , r dF(x) 
mF[T) J x ^^^^ ^r>^j,, 



xdF{x) 



F{t) > and we can write 
or equivalently 

mF{T)E[X\X <t]> 

However, E[X|X < t] < F,[X] = /i, so that the above inequality cannot hold unless mF{T) > 1. □ 

Using this lemma and Theorem 5.2 (a) (ii), the proof of Corollary 5.2 is immediate. 
Lemma 7.4. Assume r < fi and let be defined by 

r 

m 

Then, m(l - F{e)) < 1. 

Proof. As 1 — F{x) > for all a; < 6 and a.s 9 < b (because xdF{x) = ^ > ^ > 0), we deduce 
1 — F{9) > and we can write 

m(l-F(0))^°°x^^^=r</i, 

or equivalently 

to(1 - F{9))E[X\X >9]< n. 
However, E[Ar|X > 0] > F,[X] = fi, so that the above inequality implies m(l — F{6)) < 1. □ 

Using this lemma and Theorem 5.3 (a) (ii), the proof of Corollary 5.3 is immediate. 
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7.6 A remark on the unbounded case 

Recall that we supposed throughout this paper that the support of F was bounded (see assumption 
(A3)). However, viewing completeness, we should at least discuss the case of unbounded support. 

Certain cases are intuitive. For instance, it is clear that, with finite mean resource production 
and growing tails in F, the model has to allow in the limit for an infinite mean reproduction rate r. 
Otherwise survival would be impossible. (For infinite mean branching processes, see e.g. Grey [1977] 
and Barbour and Schuh [1979].) In principle, the general case of unbounded resource claims could 
be treated too, yet with some restriction on interpretability. The following theorem shows that the 
required conditions are likely to always become hard to decipher. 

Theorem 7.2. Let (S'n)„ be a sf-process on {D'^, X^, Rj^)n,k, where this double array ofi.i.d. triples 
of random variables satisfies assumptions (A1)-(A4) of Section 3 but where F has now a support 
[a, oo], a > 0. Suppose m > 1 and /i > 0. For all e > and K G M^, define pk ~ P{X > K) and 

aK,i ^ inf{i^(6) - F{a) : 0<a<b<K, b-a = e}. (77) 

Assume that, for each e > small enough, there exists some sequence {Kn)n such that 

n-PK,, and nif„(l - aK„,e)"' 0, (78) 

where px md aK,e o.^^ defined as in (77). Then, all the results of Theorem 5.3 are true. 



A new version of Theorem 2.3 of Bruss and Robertson [1991] is needed, in the case when the 
distribution is not bounded. However, Condition (78) about the tail of the distribution F is now 
needed. 

Lemma 7.5. Let Xi, X2, . . . be i.i.d. continuous real-valued non-negative random variables with 
distribution function F and mean /i < 00. Assume F having a support [a, 00], a > 0. For all e > 
and K Cz N, define pk o-nd aK,e as in Theorem 7.2, and assume that, for each e > small enough, 
there exists an increasing sequence {Kn)n of positive integers, such that 

npK„ — and nKn{l — aK„,e)'^ 0. (79) 

Let (s„)„ be a sequence such that Sn/n — > cr G (0, /i) as n ^ 00, and define 6 by 



/ xdF{x) = a. (80) 
J9 



Moreover, write M„ := Af(n, s„) :— M{n, {XkYj^^i, Sn) '^^ defined by (3). Then, 

X„_M„+i.„A0 and 1m„Ai-F(0). (81) 

n 

Proof. Let L — liminf„_j.oo ^n-Af„+i,n- Then, using the continuity of F and the Glivenko-Cantelli 
theorem. 



xdF{x) = liminf / xdF{x) (82) 



= liminf / xdFr,{x) (83) 



Sr. 



< liminf — = cr = / xdF{x). (84) 

n->-oo n Jg 

This inequality implies liminf„_j.oo -'^n-Mn+i.n > ()■ Similarly, limsup^^o^ X„_Af^_„ < 9. In order 
to prove the convergence of Xn-.M,^+i,n in probability, it thus suffices to show that X„_a/^+i.„ — 
Xn-M„,n converges to in probability. 
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Let e > 0. We write, for all K, 

P[Xn-M„ + l,n — Xn-M„,n > c] (85) 

<K]+ P[X„_M„ + l,n > K] (86) 

and we shall consider both terms separately. First, 

P[^n-A/„+i,„ >K]< P[X„.„ > A-] = 1 - (1 - pkT (87) 
where pk 1 — F{K). On the other hand, define 

a,,A' := inf{F(6) - F{a) : 0<a<b<K, b-a = f}. (88) 
Let r — \K/e\ , and /i, . . . , be a covering of [0, K] by intervals of lengths e. We have, for every m, 

P[$i,j : i^j, X,,Xj e Im] < (1 - ac^KT + n{l - ae,K)""'P[^ e /m] (89) 

<{n + l){l-a,^Kr~\ (90) 

Therefore, we get 

P[3j : Xj+i,„ - X,-„ > e] = P[3m : ^i, j, z ^ j, G /„,] (91) 

< r(n + 1)(1 - a.^Kr-' < ^^{n + 1)(1 - a,,A')"-^ (92) 
Inequalities (86), (87) and (92) give, for all n e No and for aU > 0, 

P[Xn-M„ + l,n-X„-M„,n>£] < 1 - (1 - PA')" + (« + 1) (1 - "e^i^ )"" ' • (93) 

It thus suffices to prove that the right-hand side converges to as n — > oo, if iiT is replaced by some 
sequence {Kn)n- Therefore, we need 

{l-piiS^l and 7iif„(l-ae^Aj"->0, (94) 

which is clearly satisfied if fipx^ — ^ and •nKn{\ — a^^K^)" — >■ 0, as was assumed in the statement 
of this result. 

We have thus obtained that Xn-M„+i.n 6* in probability. Therefore, we also deduce that 
M{n,Sn)/n = 1 — Fn{Xn-M^+i,n) 1 — F{9) in probability, thanks to the Glivenko-Cantelli 
theorem. □ 

Using this result, Lemma 7.2 is easily adapted to our case. 

Lemma 7.6. Let Xi, X2, ■ ■ ■ be i.i.d. continuous real-valued non-negative random variables with 
distribution function F and mean /i < 00. Assume that F has a support [a, 00], a > 0, and that, for 
each e > small enough, there exists some sequence {Kn)n such that 

npK„ and nKn{l - aK„,e)" 0, (95) 

where Pk and ax, e o^f^ defined as in Theorem 7.2. Further, let ($„)« be a sequence of integer-valued 
random variables with — 00 a.s. as n — > cx), and let (^'„)„ be a sequence of real-valued random 
variables with \l/„ — 00 a.s. as n 00. Suppose 5'n/$„ — ?► p a.s. with < p < /i and that 9 is the 
solution of 

rb 

xdF{x) = p. (96) 

Moreover, define M(t,s) := M(t, {Xk)l.^i, s) as defined by (3). Then, 

M($„, ^'„)/$„ ^ 1 - F{e), in probability. (97) 

The arguments of the proof of Theorem 5.2 are easily adapted to our case, now using Lemma 7.6 
instead of Lemma 7.2, and yield Theorem 7.2. However, as indicated before, the interpretation of 
the conditions in the unbounded case seems in general difficult. 
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8 Conclusion 



The first question posed in the Abstract is almost a philosophical one, and here Mathematics can 
hardly be of any help. We see the answer, as said in the Introduction (Section 1), somewhere in the 
middle of what contrasting theories have to offer, but this is no more than a personal perception. 
(Statistics could provide of course some answers in terms of a quantification of the likelihood of 
reasons.) 

Concerning the second question our answer is implicit in our model, and hence partially subjec- 
tive. We think that the desire to survive and to see future generations survive is the most important 
factor to keep any society together. The second important factor is in our view a sufficient standard 
of living for individuals. If individuals feel that the society could provide this but refuses to do so 
by sticking to a policy wasting resources by favouring too much certain groups, they are likely to 
emigrate. Of course, there are many other factors which may be of influence, as for instance worries 
that the society does not fight enough for human rights, or may not do enough for the environment, 
or others, but our model does not take such influences into consideration. 

The third answer is complete and interesting. The wf-society and the sf-society constitute together 
a "global" envelope in the following sense. For given distributions of natality, resource creation and 
resource consumption, little can be said "locally" , as we have pointed out before. What we find very 
interesting is that, despite this difficulty, one can show that in the long run, there is no escape, that 
is, asymptotically any RDBP (r„)„ will live between the boundaries {Sn)n and {Wn)n, and hence 
their critical parameters are of great interest. Note also that the parameters r and 6 involved in the 
computations of the critical values (see Subsections 5.3 and 5.4) may be seen as thresholds in the 
Lorenz curve we know from Economics. 

No other society in this model does as much for ensuring survival as the wf-society. The price 
to pay under the same fixed distribution is the most modest standard of living of individuals in this 
society. The sf-society forms the other extreme. Under the given assumptions this society does the 
most for the standard of living of the few. However, it jeopardises the prospects of survival severely, 
and definitely more than any other society. 

Moreover, if our perception of the weakest-first society as an extreme form of communism, and the 
strongest-first society as an extreme form of capitalism is seen as defendable, then we can conclude 
that mankind has already come close to testing the limits of societies. 



Authors' address: 
Universite Libre de Bruxelles 
Faculte des sciences 
Departement de Mathematiquc, CP 210 
B-1050 Brussels, Belgium 



27 



References 

S0ren Asmussen and Thomas G. Kurtz. Necessary and sufficient conditions for complete convergence 
in tlie law of large numbers. Ann. Probab., 8(1):176-182, 1980. 

Andrew D. Barbour and H.-J. Schuh. Functional normalization for the branching process with infinite 
mean. J. Appl. Probab., 16(3):513-525, 1979. 

N. H. Bingham and R. A. Doney. Asymptotic properties of supercritical branching processes. I. The 
Galton-Watson process. Advances in Appl. Probability, 6:711-731, 1974. 

F. Thomas Bruss. Branching processes with random absorbing processes. J. Appl. Probability, 15 
(l):54-64, 1978. 

F. Thomas Bruss. A counterpart of the Borel-CanteUi lemma. J. Appl. Probab., 17(4):1094-1101, 
1980. 

F. Thomas Bruss and James B. Robertson. "Wald's lemma" for sums of order statistics of i.i.d. 
random variables. Adv. in Appl. Probab., 23(3):612-623, 1991. 

E. G. Coffman, Jr., L. Flatto, and R. R. Weber. Optimal selection of stochastic intervals under a 
sum constraint. Adv. in Appl. Probab., 19(2):454-473, 1987. 

Harry Cohn. On the asymptotic patterns of supercritical branching processes in varying environ- 
ments. Ann. Appl. Probab., 6(3):896-902, 1996. 

Harry Cohn and Fima Klebaner. Geometric rate of growth in Markov chains with applications to 
population-size-dependent models with dependent offspring. Stochastic Anal. Appl, 4(3):283-307, 
1986. 

R. M. Corless, G. H. Gonnet, D. E. G. Hare, D. J. Jeffrey, and D. E. Knuth. On the Lambert W 
function. Adv. Comput. Math., 5(4):329-359, 1996. 

M. Gonzalez, M. Molina, and I. Del Puerto. On the class of controlled branching processes with 
random control functions. J. Appl. Probab., 39(4):804-815, 2002. 

Miguel Gonzalez, Manuel Molina, and Ines del Puerto. On L^-convergence of controlled branching 
processes with random control function. Bernoulli, 11(1):37 46, 2005. 

D. R. Grey. Almost sure convergence in Markov branching processes with infinite mean. J. Appl. 
Probability, 14(4):702-716, 1977. 

Patsy Haccou, Peter Jagers, and Vladimir A. Vatutin. Branching processes: variation, growth, and 
extinction of populations. Cambridge Studies in Adaptive Dynamics. Cambridge University Press, 
Cambridge, 2007. 

Sophie Hautphcnnc. Extinction probabilities of supercritical decomposable branching processes. J. 
Appl. Prob., 49(3);639 651, 2012. 

P. L. Hsu and Herbert Robbins. Complete convergence and the law of large numbers. Proc. Nat. 
Acad. Sci. U. S. A., 33:25-31, 1947. 

Peter Jagers. Branching processes with biological applications. Wiley-Interscience [John Wiley & 
Sons], London, 1975. Wiley Series in Probability and Mathematical Statistics — Applied Probability 

and Statistics. 

F. C. Klebaner. A limit theorem for population-size-dependent branching processes. J. Appl. Probab., 

22(l):48-57, 1985. 



28 



Karl Pearson. Tables of the incomplete beta-function. Cambridge University Press, Cambridge, 
England, second edition, 1968. 



H.-J. Schuh. A condition for the extinction of a branching process with an absorbing lower barrier. 
J. Math. Biol, 3(3-4) :271-287, 1976. 

B. A. Scvast'janov and A. M. Zubkov. Controlled branching processes. Teor. Verojatnost. i Prime- 
nen., 19:15-25, 1974. 

Kuang Xu and Shic Mannor. Rate-optimal control for resource-constrained branching processes. 

arxiv:1203.1072vl, 2012. 

Andrei Y. Yakovlev and Nikolai M. Yanev. Relative frequencies in multitype branching processes. 
Ann. Appl. Probab., 19(1):1-14, 2009. 

Nikolai M. Yanev. Conditions for degeneracy of (^-branching processes with random (p. Theory Prob. 
Applic, 20 (2):421-428, 1976. 



29 



